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ABSTRACT 

Accuracy limitations to two way range rate doi^ler tracking 
due to master (reference) oscillator frequency instabilities are 
discussed. Theory is developed to treat both the effects of 
random and non-random oscillator instabilities. The non- 
random instabilities treated are drift, environmental effects, 
and coherent phase modulation. The effects of random insta- 
bilities on range rate accuracy are shown to be describable in 
terms of cr^ (2, T, t). For the typical noise processes encoun- 
tered in precision oscillators, range rate error is related to 
the more familiar cr^ (t) and (f). Three examples are dis- 
cussed to show how to determine range rate error from given 
Oy (t) or £ (f) curves, and approximations are developed to 
simplify the treatment of complex systems. An error analysis 
of range determined from range rate data is also given. 


INTRODUCTION 

Two way range rate doppler tracking (TWDT) is a means of measuring range 
rate by observing the dcppler shift in a radio signal coherently transponded 
from a satellite. Recently proposed applications of TWDT to Earth and ocean 
physics have range rate accuracy requirements of 0.003 to 0.005 cm/ sec. ^ 
These stringent accuracy requirements impose constraints on TWDT systems 
which make careful analysis and minimization Oi system errors imperative. 

One source of error in TWDT systems is the instability of the master oscillator 
(reference oscillator). Othert have analyzed the effects of this error sourcef'^’® 
but have either not treated the problem in sufficient detail., made some errors 
in their analysis, or not treated the problem directly as it applies to TWDT 
systems for satellites. This paper will attempt to correct these defects. The 
paper will not discuss inr lvidual TWDT systems, but will take a generalized 
approach. Existing systems and requirements, however, will be kept in mind to 
ensure the utility of the paper’s results. 


1 


STATISTICS OF RANGE RATE ERROR 


A generalized schematic of a TWDT system is shown in Figure 1. hi a TWDT 
system, a reference frequency, f g, provided by tiie master oscillator is trans- 
mitted to the satellite and transponded back to the receiver. The doj^iler 
i^ifted frequency, f o + ^d* mixed with f q to produce f,}. After passing 

throu^ a band pass filter, f ^ is measured in a frequency or period counter. In 
actual systems, f q is biased by an arbitrary frequency^ so that the doi^er 
signal can be averaged over any time, T. The master oscUlator also sun;>lie8 
the time base for the coimter. The accuracy requirements for this are not 
stringent, and will not be treated here. 

The counter used to measure f^ can fit into four possible categories.^ The 
first two categories deal with the way individual measurements are taken. The 



Fig. 1— Two way doppler tracking 
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coxinter can measure the total number of periods of f^ for a fixed elapsed time. 
T, or it can measure the total elapsed time for a fixed number of periods. In 
this p^r, we shall consider T given; the results can be related to the fixed 
period case throu^ an estimate of the behavior of fi,. The other two rat^ories 
deal with the way in which successive measurements are treated. The counter 
is considered destructive or non-destructive depending on whether it is cleared 
or not cleared after each measurement. In the destructive case, the counter 
acts as a conventional period or frequency counter. In the non-destructive 
case, the counter acts as a real time clock. The effects of destructive versus 
non-destructive data taking will be discussed later. 

Nominal’y, when f q is subtracted from ail that remains is fp. The 

tracking signal, however, takes a finite time, r, to propagate from the trans- 
mitter to the receiver. If the frequency of the master oscillator shifts, in time 
T , by S fg , the counter will measure f p + 5 fg as the doj^ler shift. Since there 
is no way to distinguish between fp and S f g in the counter, b fg will introduce a 
range rate error. 

No matter how frequencies are coherently changed as signals propagate through 
the TWDT system, the doppler shift will be given by®: 

2 V 


where v is the satellite's range rate and vhere c is the velocity of light. Using 
(1), it can be shown that the range rate error from master oscillator instability 
for averaging time, T, is given by; 


6 V ( t , T, T ) 


c 

Yt 


■ ^t+T+T 

I y(t')dt' - 

-Jt +T 



( 2 ) 


y (t) is the normalized master oscillator frequency change: 

f 0 ( ^ f 0 

y(t) = -f 


and f^ is the nominal oscillator frequency. For r < T, this can be rewritten as 
(see Figure 2): 
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Fig. 2— Response functions for range rate error 


(t < T) 



c r 

pt+T+T 

M+T 1 

Sv(t,T, T) = 

y(t')dt' - 
Jt+T 

|y(f)d,'] 


If the master oscillator instabilities are random^ range rate errors can be de- 
scribed in terms of the variance of S v(t, T, t); 

<7v(T. T) = {iv(t, T, r))^^. 


< A> denotes the average of A over time, t. Using equation 3, for r < T, the 
variance becomes; 


(T < T) 



V 







where tp t,) is the two sample Allan variance of y(t) for averaging time 
tj and dead time t^ - tj given by’ ’® : 


(2. tj, tj) 


1 

2 



t + t2 

y(t')dt' 




For r > T, the variance is: 


(r > T) 

^ C7^2(2.t. T) (5) 


Notice that in the case of (5), the definitions of t and T are reversed from 
those of references 7 and 8. 

In general, <2, T, r) is not a statistic which is given in oscillator specifica- 
tions. In the time domain, generally, cr ^ iT ) = cTy*( 2 , t, t) is what is specified. 
One can relate these two statistics by’ • *: 

a^2(2,T,r) - B2(r,M)ay2(T) (6) 

where r = T/r, and where fj. characterizes the noise process involved. 

Another typical way in which an oscillator is specified is by i' (f), the power 
density in one phase modulation side band divided by the total oscillator power.® 
1’ (f) can be related to T, t) by*: 

^2(2, T,r) = - fdf £(f)|hg(f)|^ H(f) (7) 

(7T r)2 Jo 


where h^ (f - fj)) is the frequency response function of the post mixer filter of 
Figure 1, and where: 


H(f) = sin^ (tt f t) sin2 (77 f T) 
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TYPICAL OSCILLATOR INSTABIUTIES 


Master oscillator frequency instabilities can be divided into two classes, 
random and non-random effects. We shall first consider two non-random 
effects, frequency drift and coherent phase modulation. 

An oscillator's frequency drift is defined by: 

y(t) = Dt 


Using this and equations 2 and 3, one can she ^ that for both r < T and r > T: 

c cD 

8v = T ( 8 ) 

If we use 5v = 0.001 cm/sec as a target error and r = 0.3 sec as a typical 
worst case delay, we obtain for D a maximum allowable value of: 

D ^ 1.92 X 10-8 /day 


Clearly this is easily met with most good crystal oscillators. 

Equation 8 can also be used to determine environmental constraints if the fre- 
quency change per unit environmental change is known. If, for example, an 
oscillator changes its fractional frequency by 1 x 10~^ ^ per ®C, and S v < 

0.001 cm/sec is required, one can write: 


D ^ 1 X lO-^VT X K ( °C/sec) 


Using equation 8, we obtain: 


K < 0.22 °C/sec 


How coherent phase modulation of the master oscillator effects a TWDT system 
depends not only on the modulation amplitude and frequency, but also on the 
precise phase relationships between the modulation and the data taking process. 
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We can, however, obtain approxlnuite results by ignoring phase, and describing 
the modulation by: 


where f ^ is the modulation frequency, and where 1/2 is the fractional power 
in one modulation sideband. Using this in (7), we obtain: 

2Pii 


Since H (f) £ 1, we obtain the inequality: 


2P 


( 7 ^ fo ■r)’ 


Substituting this in (4) and (5), gives: 


(r < T) 


2 < 


r 


(») 


and: 


{r > T) 


2 < 


(77 fo T)‘ 


hf-(fM)r Pm 


( 10 ) 


As an example, let cr^ = 0.001 cm/sec, r - 0.3 sec, T = 5 sec, and = 
5 M Hz. We obtain: 


'hf(fM) 1^ Pm - " 112 db 
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as the estimated upper limit of filtered phase modulation that can be 
tolerated. 

Random processes typically encountered in precision oscillators can be de- 
scribed as weighted and filtered sums of power spectral densities of y (t) given 
by’-®: 

h_ ' 

S^(f) . ^ 

f 

with: 

a = 2 , 1 . 0 , - 1,-2 


Sy (f) can be related to £(f) by®: 

£(0 = 4(1^) (1®) 

Using (10), (4), (5), and previously derived equations relating cr^ (r) and 
cr^^(2, T, r ) for each noise process,® one can relate £ (f), cr^ (r) and cr^ (T, r). 
Chart 1 gives these relationships. The chart assumes that the post-mixer band 
pass filter is infinitely sharp with a pass band of fjj - f^ to f q + fj^ , that 
2 77^ fi^ T 1, and that 2 77 f^ T » 1. For completeness, drift is included even 

though it is not a random process. To simplify formulae in the chart, the fol- 
lowing are used: 

_ T 
r - — 

r 


X - J y d t 

C l for r - 1 
^ loforrT^O 


phase 
2 77 
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Chart 1— Range Rate Error for Various Noise Processes 







Notice that for all processes but random walk y and drift, when r < T, Is 
proportioikal to By statistically averaging N successively taken counter 
readings, one improves the variance of the average by or This means 
that as long as random walk y or drift processes a.'e not involved, smaller 
errors are obtained by going to longer averaging times or taking non-destruc- 
tive data rather than by statistically averaging dr.ta. Even when taking non- 
destructive data, to obtain the smallest cr^ for 'x given T, one should only use 
data every T seconds apart, and not least squrxes fit the data from smaller 
time intervals. 

Notice also that for white x noise, is (roportional to f^. This allows one to 
reduce the white noise contribution to by reducing f Flicker x also de- 
pends on f^, but in such a slowly varying manner that it can essentially be con- 
sidered fixed. 

APPUCATION TO TYPICAL REFERENCE OSCILLATORS 

In this section, cr^(T, t) shall be derive^I for three reference oscillators, an 
Oscilloquartz crystal oscillator, a Hewlett-Packard cesium staixlard, and a 
NASA hydrogen maser, to demonstrate techniques for using the theory of the 
previous sections. The results for the two commercial oscillators used are 
based on manufacturers specifications. Their use is for the {mrpose of example 
only, and does not constitute an endorsement of these products or a confirmation 
of their specifications. 

The first oscillator we will consider is the Oscilloquartz B-5400 crystal oscil- 
lator.^® Its ii’(f) spectrum is given in Figure 3 and its <^yir) curve is given in 
Figure 4. Deriving for this device is very simple; since its instabilities are 
just the sum of simple processes (white y, flicker x, flicker y, drift), is 
just tlie sum of the results taken from Chart 1 for each process. Figure 5 
shows o^(T, r) for t = 0.3 s. 

The second oscillator to be considered is the Hewlett-Packard 5061A High Per- 
formance Cesium Beam Standard. ^ (f) for this device is shown in Figure 6. 
Notice that, in this case, we don't have just the sum of simple processes. 
Equation 7 could be used to generate from £ (f), but this tedious method can 
be avoided; we can quickly obtain an approximate curve for <r^ (T, r) by using 
the limiting properties of H (f) in relation to £ (f). 

For our approximation, we shall rely on the fact that £ (f) can be broken in two 
general categories: a short term £ (f) which changes slowly with f, and a long 
term £{> (f) which blows up as f goes to zero and makes a negligible contribution 
for large f. Using this, we can write; 
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Fig. 4— cTy (t) for Oscilloquiartz B-5400 



(cm/s) 
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cr 2 (2 T, r) ~ cr 2 + cr .2 

y \ ' y ys yt 

and use (4) and (5) to obtain approximations for 

Since H(f) is a periodic fimction, for a slowly changing (f), we can replace 
H(f) £-(f) in the integral of equation 7 by < H(f) > £^(f). < H(f) > denotes the 
average of H (f). The short term part of equation 7 becomes: 


r) 


^ ■^1 ^ r . 1 




Jjh£(f) |^ £ (f) df 
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Fig. 6— 5^(f) for Hewlett-Packard High Performance Cesium Standard 


The long term part of cTy^(2, T, t) is determined by H (f) (f) in a r^ion near 

f = 0. For long term noise processes, T, t) is determined by H(f) ££(f) 
in the regions of f given in the following chart (r < T); 


noise nrocess 


region of contribution 


flicker y 


random walk y 




For a complex (j[)» which jmlae process contributeB is determined by the 
chart and the appropriate values of T and r. For the Hewlett-Packard cesium 
standard* £ (f) changes from flicker y to white y noise at f = 0.1 Hz. Ihis 
means that, for -r - 0.3 s, is determined by the flicker y part of £(f) for T 
up to 5 seconds, and by the white y part of £ (f) for T greater than 5 seconds. 
Because white y noise is determined by £ (f) for f up to 1/2 r = 1.67 Hz, and the 
white y noise is truncated at f = 0.1 Hz, the full value of cr^ from Chart 1 can- 
not be used to determine the vdiite y contribution for the cesium standard, but 
must be multiplied by a reduction factor, p. An approximate method for calcu- 
lating p is shown in Figure 7. Figure 8 shows (T, 0.3 s) for the cesium 
standard using the approximation techniques just derived. 
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The last device we will consider is NASA's latest generation hydrogen maser. 
The ^y("^) curve is shown in Figure 9. The device uses an Oscilloquartz B-5400 
as a local oscillator which it phase locks to the maser with a 1 Hz loop. Using 
this fact and Figure 9, by indirectly using (11), we can calculate the white x 
noise contribution to a-J with a 1 KHz band width as the sum of the oscilloquartz 
contribution with a 1 KHz bandwidth and the maser's contribution with a 1 Hz 
bandwidth. In this case, the other noise processes make a negligible contribu- 
tion. cr^ (T, 0.3 sec) is shown in Figure 10. Even though is determined by 
white X noise, cr^ cannot be appreciably reduced by reducing f^^; is deter- 
mined principally by the maser white x noise in the 1 Hz loop bandwidth. 

RANGE FROM RANGE RATE 

Integrating Iq yields the integral of range rate: range. Since the averaging 
process involves intpgj ation, not dividing the information stored in the counter 
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register by T tiirns range rate data into range data. This means that the error 
in taking range data produced by the master oscillator is just: 


S s(t, T, t) - T Sv(t, T, t) 


Similarly, this also means for random processes; 

or/(T, t) = 


For the three sample oscillators, Figure 11 shows ^^(T, 0.3 s). Not'ce again 
that taking non-destructive data instead of destructive data will yield oetter re- 
sults in all cases but that of the crystal oscillator for T greater than 400 sec- 
onds. 
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